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ABSTRACT 

Fluctuation properties of squeezed photon beams gener- 
ated in three wave mixing processes such as second har- 
monic generation, degenerate and nondegenerate para- 
metric oscillations, and homodyne detection are studied 
in terms of photon sequences recorded by a photode- 
tector. 


Photon number fluctuations and photon number 
correlations are fundamental properties of a light beam. 
These properties are different for different light sources 
and can be used to characterize photon beams. In this 
short communication we discuss our work on statisti- 
cal properties of squeezed photon beams generated in 
three wave interaction processes in terms of counting 
and waiting time distributions. We summarize some of 
the interesting results obtained for these systems. Pro- 
cesses that we consider here include second harmonic 
generation, and degenerate and nondegenerate paramet- 
ric down conversion (DPO and NDPO). Squeezed state 
of light have been realized in these systems experi- 
mentally (Rcf.l). Homodyne statistics when squeezed 
light produced by the DPO is mixed with coherent light 
from a local oscillator arc also discussed. A dynamical 
model for these beams is used and photon sequences 
recorded by a photodetcctor are calculated. 

We use positive-P representation (Ref.2) to map 
quantum mechanical equations of motion for the anni- 
hilation and creation operators onto a set of C-number 
stochastic equations for the complex field amplitudes. 
Using simple transformation of field variables it can 
be shown that the field produced in these processes 
can be described in terms of independent real Gaussian 
stochastic processes (Ref.3-4). 

We use a generating function technique to ob- 
tain the statistics of the photons emitted by these light 
sources. The generating function G\.s, t, T) for the pho- 
ton statistics measured by a detector with a parameter 
.s is given by (Rcf.5) 


G(s,t,T) 



I(t')dt' 

>• 


(i) 


Here ?/ is detector efficiency and /(/) is photon flux 
emitted by the source. Generating function G(\,t,T) 


is simply the probability of detecting no photon in the 
time interval Mo / + T. In order to obtain generating 
function we express I(t) in terms of the c-number field 
variables. Statistical averaging is performed by mak- 
ing Karhunen Lo£ve expansion of the field variables in 
terms of a set of orthogonal functions. Following the 
method developed in our earlier investigations (Ref.3) 
we derive an analytical expression for the generating 
function G(s,t,T) for the photon statistics. From this 
generating function various statistical quantities such as 
factorial moments, photon counting and waiting time 
distributions can be obtained. 

The photon counting distribution p(mJ,T) is 
the probability of counting m photons in the time in- 
terval t to / +T. It can be obtained from the generating 
function by using the relationship 


]AmJ,T) 


(- l) m 


d m 

ds ] 




j=i 


(2) 


The waiting time distribution w(t,T) is the probability 
density for two successive photoelectrons to be sepa- 
rated by the time interval 7 and it is given by 

d? 

w(t, T = <' — <) = -t). (3) 

dtdt f 

In the stationary regime these quantities are indepen- 
dent of initial time t. Here we only summarize photon 
statistics only in the steady state regime. The field from 
the DPO can be expressed in terms of two independent 
Gaussian random variables with mean zero and differ- 
ent variances 

< UiiDUjU + T) >= &ij c~ x,T . (4) 


Here k is the mode coupling constant, and e is the 
dimensionless amplitude of the pump beam incident on 
the cavity. The decay constants Ai and A 2 are given 
by 

Ai = (7 - |« f |) , A 2 = (7 + \k(\) . (5) 

Here (I/27) is the cavity lifetime. Below' threshold 
Aj and A 2 are always positive. Using the properties 
of field variables ui and u 2 we obtain the generating 
function for the DPO as (Ref.3) 


G(s, T) = Q\ (s, T)Q?(s, T ) , (6) 
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where 


Qi{s,t,T) 


[cosh(ijT) + fi(l) sinhfoT)] 


(7) 


with 


and 




( 8 ) 


z\ = + 2si]jkc, 

Z 2 = ^2~ 2st]yK(. 


(9) 


Mean photon number inside the cavity is given by 




( 10 ) 


For the DPO £ is equal to half. Once the generating 
function is known photon counting and waiting time 
distributions are obtained from Eqs. (2) and (3), For 
small mean photon number and short counting time 
p(m,T) decreases monotonically. 



Figure la. Photon counting probability distribution for 
the DPO for R = 0.95, unit efficiency, and counting 
time interval 2y T = 20 . 

For long counting times p(m,T) shows sharp 
even-odd oscillations implying that the probability of 
detecting odd number of photons is much smaller than 
the probability of detecting even number of photons. 
As the mean photon number h is increased these even- 
odd oscillations become smaller. Figure (la) shows 
p(rn,T) for R = k c jy - 0.95 near threshold corre- 
sponding to n = 4.63. These curves are meaningful 
only for integer values of m. We see that near thresh- 
old even-odd oscillations become less pronounced and 


p(rn,T) develops a long tail. We have also studied 
photon statistics for the DPO in the transient regime, 
that is, during its evolution from vacuum state to the 
steady state (Ref.4). For small transient time even-odd 
oscillations in photon counting distributions are even 
sharper than the even-odd oscillations in the stationary 
regime. 

With the degenerate modes of a parametric os- 
cillator, nondegenerate modes are also present. We 
consider the nondegenerale modes of parametric os- 
cillator for which two nondegenerate photons have the 
same frequencies. These fields can be expressed in 
terms of four real random Gaussian variables. Here 
we discuss two cases, one in which amplitudes of the 
two nondegenerate modes are homodyned and second 
in which intensities of the two modes are added to- 
gether. For a given pump strength the NDPO the mean 
photon number is much smaller than the mean photon 
number for the DPO. Mathematical expression for the 
generating function for the first case is similar to the 
generating function for the DPO. Thus p(m,T) for the 
NDPO also shows even-odd oscillations. However, for 
the same pump strength, even-odd oscillations in the 
NDPO are sharper then the oscillations for the DPO. 
They are centered towards smaller values of m. The 
difference between the DPO and the NDPO lies mainly 
in the value of the mean photon number. 


CM 



Figure lb. Photon counting probability distribution for 
the NDPO for R = 0.95, efficiency t] - 1 , and counting 
time interval 27 !’ = 20 . 

The second case that we consider is when inten- 
sities of the two modes are added together. In this case 
the power £ appearing in Eq. (7) is one (Ref. 6 ). In this 
case the expression for Q x ( $ , T) is similar to the gener- 
ating function for thermal light (Ref. 5). The counting 
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distribution for NDPO, however, is very different from 
that for the thermal light. It shows even-odd oscilla- 
tions as a function of m whereas such oscillations are 
not seen for thermal light. Figure (lb) shows even- 
odd oscillations in p(rn, t) for R = k c/y = .95 for the 
second case. 

Next we discuss photon statistics of the funda- 
mental beam from an intracavity second harmonic gen- 
eration (SHG). Field for this system can be expressed 
in terms of two real Gaussian random variables and a 
coherent component (Ref.7). These results are obtained 
by linearizing the field amplitude equations around the 
deterministic steady state values. The generating func- 
tion for the SHG can be written as (Ref.6) 

G(s,T) = Qi(s,T)Q 2 (s,T)e- J(, T) . (10) 


Here 

/(s, T) = -2 ST}hn 0 T 
2sr]hv(zi ) 


H( l+ 2 

z\ V (2 + AjT) 


+ ( 2 + A iO(i + ^) 


(ID 


with 


v(z l) = 


1 + 


( 2 ty ) tanh(ziT/2) 
1 + ^ tanh(ziT /2) 


( 12 ) 


Here Qi(s,T) and Q 2 (s,T) are given by equation (7 ) 
with / = 0.5 and 


A i = (1 + 3n) , Ao = (1 + n) , (12) 


For the SHG z i and z 2 are given by 

z\ = \\ — 2si]fi , 
z\ = A 2 + 2srjh. 


(13) 


Here h and n 0 are the average and threshold pho- 
ton numbers, respectively. From this generating func- 
tion various statistical quantities of interest can be cal- 
culated. Photon sequences in the SHG can be an- 
tibunched. Although the antibunching effect is very 
small riding on an intense coherent background it is 
clearly reflected in the behavior of the waiting time 
distribution for the SHG. 

Using similar techniques we also obtain the gen- 
erating function when light from the DPO is homo- 
dyned with coherent light from a local oscillator. De- 
pending upon the relative phase between coherent light 
from the local oscillator and squeezed light from the 
DPO we can see sub-Poissonian or super-Poissonian 


statistics for the homodyned photon beam. Figure (2) 
shows waiting time distribution when the relative phase 
is 0° and 90°. We see bunched light when coherent 
component is added to the unsqueezed component and 
antibunched light when coherent component is added 
to the squeezed component. 
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Figure 2. Waiting time distribution for DPO mean pho- 
ton number ft = 0.2 and local oscillator mean photon 
number n 0 = 2. Dashed curve is for coherent light 
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